global soal1 soal2
global numsoal
global x0 xn h n y01 y02 numy0
global fungsi1 fungsi2
global numfungsi
global metodestring
global hasil1 hasil2 
global numhasil

x = [1900 2000];
n = 10;
h = (x(2) - x(1)) / n;
y0 = 76.4; % * 10^6; % 76.4 juta jiwa
y_absen = @(x,y)(0.02*y-0.00004*y^2); % P itu x atau y?
orde = [ 3 4 ];

% Orde 3
[n result1] = Runge_Kutta( x, h, y0, y_absen, orde(1));

% Orde 4
[n result2] = Runge_Kutta( x, h, y0, y_absen, orde(2));

result = [result1 result2];

clear result1 result2;

%------------------------------------------------%
% Pembuatan report
%------------------------------------------------%
% Soal
soal1 = 'P'' = 0.02*P-0.00004*P^2';
numsoal = 1;
% Variabel
x0 = x(1); xn = x(2);
% h = h (variabel dan globalnya sama)
% n = n (variabel dan globalnya sama)
y01 = y0;
numy0 = 1;
% Fungsi untuk perhitungan
fungsi1 = func2str(y_absen);
numfungsi = 1;
% Metode
metodestring = 'Runge-Kutta Orde 3';
% Hasil
hasil1 = result(1);
numhasil = 1;
% Generate Report
report('DifferentialEquation');
